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We present a theory for the elasticity of cross-linked stiff polymer networks. Stiff polymers, unlike 
their flexible counterparts, are highly anisotropic elastic objects. Similar to mechanical beams 
stiff polymers easily deform in bending, while they are much stiffer with respect to tensile forces 
("stretching"). Previous approaches have based network elasticity on the central- force stretching 
mode, in a manner similar to classical rubber elasticity for flexible polymers. In contrast, our 
theory properly accounts for the soft bending response inherent to any stiff polymer network. A 
self-consistent effective medium approach is used to calculate the macroscopic elastic moduli starting 
from a microscopic characterization of the deformation field in terms of "floppy modes" - low-energy 
bending excitations that retain a high degree of non-affinity. The length-scale characterizing the 
emergent non- affinity is given by the "fiber length" If, defined as the scale over which the polymers 
remain straight. The calculated scaling properties for the shear modulus are in excellent agreement 
with the results of recent simulations obtained in two-dimensional model networks. Furthermore, 
our theory can be applied to rationalize bulk rheological data in reconstituted actin networks. 

PACS numbers: 62.25.+g, 87.16.Ka, 81.05.Lg 



I. INTRODUCTION 

The elasticity of flexible polymer gels is successfully 
described by the theory of rubber elasticity [lj. It as- 
cribes the resistance to deformation to a reduction of 
conformational entropy induced by a changing end-to- 
end distance of individual polymer strands. In the classic 
approach, developed by Kuhn and others [2[ , the magni- 
tude of the deformation of a single constituent polymer 
is usually assumed to derive from the macroscopically in- 
duced strain in an affine way. With this assumption the 
network problem is reduced to calculating the response 
of a single chain. In this sense affine deformations repre- 
sent a mean-field assumption that neglects spatial corre- 
lations and therefore the coupling between the network 
structure ( "architecture" ) and the mechanical properties 
of its constituents. 

In recent years a different class of cross-linked net- 
works made of semiflexible or stiff polymers have gained 
widespread interest. Their importance for biological sys- 
tems as the cytoskeleton or extra-cellular matrix makes 
understanding their properties highly rewarding [H . Out 
of the variety of biological stiff polymers, F- actin has 
emerged as a model system, which allows precise in vitro 
rheological measurements, for example in determining 
the (complex) frequency-dependent shear modulus G(o>) 
and in particular its elastic component, the plateau mod- 
ulus Go at intermediate frequencies. In these experi- 
ments various types of cross-linking proteins are being 
used 0, H, H, 0] and the influence of the degree of cross- 
linking on the elastic modulus is investigated. 

Stiff polymers, unlike their flexible counterparts, are 
highly anisotropic in their elastic response and may be 
characterized in terms of two qualitatively different defor- 
mation modes (see Fig. [1} [1, Q . The linear response to 



longitudinal forces acting parallel to the contour (stretch- 
ing/compression), is due to the presence of thermally ex- 
cited undulations similar to the (isotropic) stiffness of 
flexible polymers. The resulting effective spring constant 
of a stiff polymer of contour length l s , kn ~ de- 
pends on the temperature-dependent persistence length 
l p ~ T -1 , which indicates the entropic origin. On the 
other hand, the resistance of the polymer to transverse 
forces (bending) is predominantly an energetic effect, 
leading to an increase in energy rather than to a decrease 
in entropy. Subsequently, the corresponding spring con- 
stant k± ~ lj 3 is independent of temperature. 
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FIG. 1: Sketch of a stiff polymer network with filaments which 
are straight on a scale If and where the distances between 
crosslinks ("polymer segments") on a given filament is de- 
noted by l 3 - The response of the constituent stiff polymers 
to external forces is anisotropic with spring constants ku and 
k±, characterizing their resistance to stretching and bending 
deformations, respectively. 

The presence of two elementary deformation modes 
complicates, but also enriches, the theoretical analysis of 
stiff polymer networks since it is not obvious which of the 
modes, or combination thereof, will dominate the macro- 
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scopic elastic response [1CJ, |ll| . Leaving aside these dif- 
ficulties, recent approaches |9l [l2l. [l3j have still adopted 
straightforward extensions of rubber elasticity to stiff 
polymer networks by assuming affine deformations to be 
present down to the scale of the individual polymer seg- 
ment - the part of a polymer filament that connects two 
neighbouring cross-links (see FigfT]). In these models no 
bending deformations are present, leaving the stretching 
mode as the only possible source of elasticity. 

In many systems of interest, however, the large value of 
the persistence length l p /l ;§> 1 calls this affine approach 
into question. This separation of length scales implies 
that the bending mode is in fact soft as compared to the 
stretching mode, since k^/k± ~ l p /l ^> 1. One would 
therefore expect the elastic energy to be dominated by 
low-energy bending deformations instead of highly ex- 
pensive stretching modes [lij]. Interestingly, recent sim- 
ulations on random fibrous networks have shown that this 
is not always the case [ill HH, EH EH E3|- There, it was 
found that in networks with infinitely long filaments (for 
fixed density) the soft bending mode is suppressed and 
the elastic modulus is entirely given in terms of the stiffcr 
stretching mode, similar to springs connected in paral- 
lel. In contrast, the same simulations performed in the 
more realistic situation of finite filament length have in- 
deed identified an elastic regime that is dominated by soft 
bending deformations. The filament length thus strongly 
influences the elastic properties and is crucial for under- 
standing the observed behaviour. The affine theory, by 
working on the smaller scale of the polymer segments, is 
incapable of accounting for these effects. 

In the present article we expand on our recent publi- 
cation [18| to develop an elastic theory that works on the 
scale of the whole polymer filament. The theory naturally 
explains the presence of a bending dominated regime as 
well as its suppression with increasing filament length. It 
is expected to be applicable to a broad class of filamen- 
tous networks with a soft bending mode. 

Similar to the classic theory of rubber elasticity it 
assumes that the cross-links adjust to the macroscop- 
ically applied strain without showing thermal fluctua- 
tions. In contrast to rubber elasticity, however, the cross- 
link movemens are chosen such that the polymer end-to- 
end distances are kept unchanged. This automatically 
avoids energetically highly expensive stretching deforma- 
tions and results in elastic moduli that derive from the 
soft bending mode only. 

In the following, we assume that stiff polymers, char- 
acterized by kn/k± ^> 1, effectively behave as if they 
were strictly inextensible bars, i.e. having an infinite 
stretching stiffness fcii — -> oo. Subsequently, we will con- 
struct sets of "admissible" cross-link displacements that 
respect this inextensibility and thus retain a highly non- 
affine character. These displacement modes are referred 
to as "floppy modes" [ll|, highlighting the fact that in 
an equivalent network of central force springs they would 
carry no energy. Here, the finite bending stiffness of 
the polymers associates an elastic energy to each mode, 



which we use to calculate the macroscopic elastic con- 
stants of the network. 

Section |TT] will be concerned with the analysis of net- 
works in the limit of diverging stretching stiffness fc|| — > 
oo, which allows us to treat stiff polymers as inextensible 
bars. We will introduce the concept of the floppy modes 
and give an explicit construction valid for a broad class 
of network architectures. 

In Section [III] we discuss the energy involved with ex- 
citing floppy modes in networks of stiff polymers, char- 
acterized by a finite, but soft, bending stiffness. Specif- 
ically, we will develop a theory that allows to calculate 
the network elastic constants in a self-consistent manner. 

Section IIVI is devoted to the specific architecture of 
random fibrous networks in a planar geometry (two di- 
mensonal), where we check our ideas against simulations. 

II. FLOPPY MODES 

Here, we are concerned with some general properties 
of networks of inextensible bars, so called frameworks. 
While the bars are assumed to be perfectly rigid, they 
are allowed to freely rotate at the cross-links ("vertices"). 
In effect, both the stretching and the bending mode are 
eliminated, which leaves us with a purely geometric prob- 
lcm. By applying methods from rigidity theory [1 91 ] we 
will find that polymer networks when viewed as frame- 
works are not rigid and possess zero-energy deformation 
modes ("floppy modes"), for which we will give an ex- 
plicit geometric construction. These modes, which may 
be viewed as the analog of the zero-energy shear modes 
of regular square lattices, characterize the deformation 
field of the network under external strain. By account- 
ing for the finite bending stiffness of the polymers, they 
are used to calculate the elastic energy stored in polymer 
networks and thus the elastic moduli. 



A. Maxwell counting 

It has first been realized by Maxwell [2(| that a frame- 
work, consisting of v vertices and b bars, can undergo 
a transition from a floppy to a rigid state by increas- 
ing the coordination number z. Assuming that each bar 
represents an independent constraint for the total of dv 
degrees of freedom in d spatial dimensions, Maxwell de- 
rived the condition b — dv = determining the rigidity 
transition. As the number of vertices can be rewritten 
in terms of the coordination number as v — 2b/ z, this 
immediately yields a critical coordination of z c = 2d. 
According to this simple Maxwell counting rule, frame- 
works are rigid, whenever their vertices have more than 
z c neighbors, while they will be floppy and allow for in- 
ternal rearrangements otherwise. 

With regard to stiff polymer networks this transition 
may be used to set up a classification where the elas- 
tic energy is dominated by either bending or stretching 
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modes. While for z < z c bending modes can stabilize 
the otherwise floppy (zero-energy) central-force network, 
they only provide minor contributions to the energy once 
z > z c . The honeycomb lattice in 2d, for example, has 
a coordination of z = 3 and is therefore bending domi- 
nated, while the triangular lattice with z = 6 is clearly 
rigid and therefore stretching dominated. Imposing a de- 
formation necessarily leads to the stretching of bonds. A 
particular case is the square lattice in two dimensions, 
which has precisely the critical coordination z = z c = 4. 
Although being floppy with respect to shear deforma- 
tions, the network may be stabilised by introducing suit- 
able boundary constraints or by adding additional bars 
along the diagonals of some of the squares. It turns out 
that in the limit of infinite system size a non-extensive 
number of diagonal bars (which scales as the square root 
of the system size) is needed to stabilise the network [2l[ . 

Maxwell-counting is only approximate, since one can 
always add redundant bars that do not constrain any de- 
grees of freedom. This effect is taken into account by the 
modified Maxwell relation b — dv = s — m (22J. In this pic- 
ture redundant bars create overconstrained regions where 
a total of s states of self-equilibrated internal stresses may 
exist. In general, a state of self-stress is defined as a set 
of bar tensions that is in static equilibrium with zero ex- 
ternal force applied. At the same time underconstrained 
regions arise that allow for to zero energy deformation 
modes, i.e. internal rearrangmcnts that can be accomo- 
dated without changing the lengths of any of the bars to 
first order in the magnitude of the imposed strain. These 
are usually referred to as mechanisms or floppy modes. 

In principle, the floppy modes of a pin-jointed struc- 
ture may be found by studying the kinematic matrix C 
which relates vertex displacements d to segment exten- 
sions e = Cd [23| . The kinematic matrix thus constitutes 
a linear relation between displacements and extensions, 
which is only true for infinitesimally small displacements. 
The entries to the matrix can then easily be identified by 
considering the extension of a single bar oriented (in two 
dimensions) at an angle <f> to the horizontal. For given 
displacements = (ui,Vi) at the two vertices i = 1,2 
the extension is found as 



(u 2 - tti ) COS (j) + (V2 -Vx)\ 



(1) 



The floppy modes then correspond to those vertex dis- 
placements that do not lead to any extensions in the bars. 
This amounts to calculating the null-space of the matrix, 
i.e. Cd = 0. 

An elementary but illustrative example of a bar/joint 
network (adopted from Ref. [24[) is the "chair" shown 
in Fig. [2^. Having b — 4 bars and v — 2 vertices, 
Maxwell's counting rule would imply that the structure is 
marginally rigid. Actually, there is also one floppy mode 
rn = 1 as well as one state of self-stress s — 1 . The former 
corresponds to the (infinitesimal) movement of the hori- 
zontal bar forming the seat, while the latter corresponds 
to a tension in the two vertical bars making the back. 

For regular systems it is sometimes possible to guess 
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FIG. 2: Illustration of the floppy modes. The "chair" in a) has 
one floppy mode (FM) corresponding to the axial movement 
of the horizontal bar, as well as one state of self-stress (S) 
located in the two vertical bars. The floppy modes of the 
honeycomb lattice may be constructed, b) from the global 
shear deformations along any of the three dashed lines (as 
well as their parallels), or, c) from localized librations. 



the modes. Consider, for the purpose of illustration, a 
honeycomb lattice in two dimensions, where a coordi- 
nation of z = 3 implies b — 2v — —v/2. There is, ac- 
cordingly, half a floppy mode per vertex. These modes 
are most easily identified with shear deformations along 
lines of symmetry (Fig. Wp)- Probing the shear response 
of the honeycomb along a given direction will cause each 
of the N layers of cells to be displaced by a small amount 
(5, which eventually has to add up to the externally im- 
posed deformation A = NS. Thus, there is "sharing" of 
the deformations between the individual cells and each 
layer contributes a small amount to fulfilling the con- 
straints imposed by the macroscopical strain field. In 
other words, the deformation field in the honeycomb lat- 
tice is affine down to the scale of the individual cell, which 
experiences deformations 6 — A/N cx l cc u proportional 
to its own size. 



Another possibility to construct the floppy modes of 
the honeycomb network is given by the librations of indi- 
vidual hexagons [25( (see Fig. [2J3). These librations are, 
in contrast to the shear displacements, localized modes 
that are confined to a single cell and its immediate sur- 
roundings. Since there is one libration per cell and each 
of the six corresponding vertices belongs to three cells, 
this also makes one mode for every two vertices. 
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B. Floppy modes of stiff polymer networks 

Proteins used to cross-link stiff polymers into networks 
often have only two heads [26| such that there can only 
be two-, three- and four-fold connected vertices. The 
average coordination number in stiff polymer networks is 
therefore z < 4, which would place it below the rigidity 
transition and render the network bending dominated. 

In contrast to the very regular structures discussed 
above, stiff polymer networks are usually highly random. 
Nevertheless, as we will see below, the floppy modes can 
be constructed quite easily on scales If over which the 
undeformed polymers can be assumed to be represented 
by straight fibers. For isolated polymers the length-scale 
If can be identified with the persistence length l p , while 
in networks the origin may be different and for example 
a consequence of the network generating process itself. 
It is the presence of the length-scale I / which renders the 
structure of stiff polymer networks qualitatively different 
from flexible polymer gels. The resulting fibrous appear- 
ance may be inspected in the figures of Refs. [13, Hg, [2i| . 
We have recently argued that the "fiber length" I / plays 
the role of the size of an effective unit cell [lj| . This has 
to be contrasted to flexible polymer gels, where the unit 
size is set by the mesh-size. In the following we use the 
word "fiber" in connection with the length I f over which 
the polymer remains straight. In later sections we will in- 
troduce a simple model system where fiber and polymer 
length are equal. 

Consider (see Fig. [3K,b) a single (primary) fiber a 
of length If, which may be part of a longer polymer. 
It is imbedded into a network of other (secondary, ter- 
tiary, . . .) fibers. Secondary fibers a* are assumed to in- 
tersect the primary fiber a at the crosslinks i — 1 . . . n c \, 
while tertiary fibers only intersect secondary fibers, e.t.c. 
The floppy-mode construction proceeds in two stages 
during which only the cross-links on the primary fiber 
are being moved. The rest of the network, in particular 
the neighbouring secondary filaments will remain static 
such that the floppy mode stays highly localized similar 
to the librations of the hexagons discussed above. In the 
first step, we perform a small axial displacement Sz a of 
the primary fiber a as a whole. The axial movement of 
the cross-links pertaining to this fiber induces a change in 
length of all neighbouring segments on the crossing sec- 
ondary fibers. In the second step, therefore, one has to 
account for the length constraints on these segments by 
introducing cross-link deflections y a ^ transverse to the 
primary fiber. It turns out that to first order in Sz a all 
segment lengths can be kept at their repose length by 
choosing 

y a ,i = Sz a cot 6> Qjl e Qjl + , (2) 

where 8 a ^ is the angle between the two crossing fibers a 
and at at crosslink i. We denote by e a ^ a unit vector 
transverse to the primary fiber lying in the plane spanned 
by the two fibers, and by y^ i an arbitrary vector per- 
pendicular to this plane. 




FIG. 3: Construction of a floppy mode starting from the 
initial geometry as drawn in a). In b), the horizontal fiber 
is moved, while the surrounding fibers remain in their orig- 
inal positions. This leads to the new cross-link positions 
(green circles) with transverse deflection y a] i — — cot 9 a< i Sz a 
(Eq. ((2J) ). The component y^,i (not drawn) is oriented per- 
pendicular to the plane of the two fibers. In c), also the 
secondary fiber is moved, such that the cross-link is now de- 
flected according to y a ^ = y a ^ + Sz ai / sin# aj i (Eq. ©). The 
solid green line represents the actual contour of the deformed 
fiber obtained by minimizing the bending energy along the 
entire fiber (see Eq. ((5JI). 



We would like to emphasize that the construction only 
works for infinitesimal Sz a , while finite displacements nec- 
essarily lead to changes in bond lengths and therefore to 
stretching of bonds. As will be explained in more de- 
tail in Section IIV Al this has dramatic consequences on 
the nonlinear elasticity of the network, leading to strong 
strain stiffening behaviour. 

The construction can be performed for any of the 
a = 1 , . . . , Nf fibers, such that precisely Nf floppy 
modes are identified this way [42 |. For the mode lo- 
calized around fiber a one may define a vector y a — 
(0 ■ ■ ■ , y a ,i, • ■ ■ , y a .n cl , 0, . . .), where the deflections of all 
crosslinks in the network are combined. It has nonzero 
components only at crosslinks belonging to fiber a. With 
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respect to the standard vector scalar product one can 
then show that the set of floppy modes {y a } is linearly 
independent, however, not orthogonal. Since a given 
crosslink i always belongs to two filaments at the same 
time, there is obviously a coupling between the two cor- 
responding modes. 

In analogy to the shear modes of the honeycomb lat- 
tice one may also construct extended floppy modes for 
the fibrous network by superposition of the different y a . 
This amounts to relaxing the constraint that a single fiber 
moves in a static environment where neighbouring fibers 
remain fixed to their initial positions. Instead, a different 
Sz a is assigned to each of the a — l,...,Nf fibers. In this 
case, we find (see Fig. |3J:) that the crosslink deflection 
y a .i of Eq. @ has to be modified by a term Sz a J s'm9 a ,i 
due to the additional movement 5z ai of the neighbouring 
filament a 2 : at crosslink i. This amounts to the overall 
deflection 

5z a . 

Ya.i — Ya.i "i" ~. 7 &a A • (3) 
Slllffa,, 

For the particular architecture of a random fiber net- 
work in two dimensions (2d), to be introduced below (see 
Section lTV)) . we have also obtained an orthornormal set of 
floppy modes. The appropriate values of y a ,i and 5z a are 
found by performing a singular value decomposition of 
the compatibility matrix C. One of the modes is visual- 
ized in Fig. 21 where the black lines indicate the floppy- 
mode displacements of the crosslinks. One remarkable 
property is the heterogeneous distribution of amplitudes 
x, which leads to polynomial tails in the probability dis- 
tribution, P(x) ~ a; -3 (see Fig. [5]). The exponent is a 
direct consequence of the random orientation of the fila- 
ments which induces a probability distribution of angles 
9 between two intersecting filaments, P{9) ~ sm.9 [43]. 
By a transformation of variables to the floppy mode de- 
flection y ~ cot 9 (and thus to x) one finds a distribution 
P{y) ~ sin 3 9(y) — > y~ 3 , where the latter limit corre- 
sponds to large y 1. 



III. ENERGY OF FLOPPY MODES 




FIG. 4: Floppy mode of a random fiber network with 225 
fibers (see Sect llVfl . A fiber has the length of one third of the 
system size. Grey lines represent the network, black lines the 
floppy mode displacements. Note, that the overall ampitude 
of the mode is arbitrary. 
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FIG. 5: Normalized probability distributions of absolute val- 
ues of floppy mode displacements, as shown in Fig. 3] The 
distributions for different fiber numbers Nf can be rescaled 
on a single master curve by changing the overall amplitude of 
the modes. 



Having constructed the floppy modes of the fibrous 
polymer network we now proceed to determine the energy 
associated with their excitation. 



A. Effective medium theory 

We calculate the response of the network to the axial 
movement Sz a of only one fiber a. This proceeds by first 
calculating the energy stored in the primary fiber, given 
that the surrounding fibers are not allowed to move. Sub- 
sequently, we will relax this constraint and also account 
for the energy induced in the secondary (tertiary, . . .) 
fibers. We will find that infinitely many levels of neigh- 



boring fibers can be included by a formal resummation 
of the energy terms in a Cayley-tree approximation. 

Quite generally, the bending energy of a polymer fiber 
with weakly undulating contour can be written as 

where y Q = y(s a ) denotes the transverse deflection at 
point s a along the backbone of the polymer a. The 
bending rigidity n is related to the peristence length by 
k — IpkgT. The actual value of the energy contained in a 
localized floppy mode of amplitude Sz a can then be found 
by minimizing Wb for the given set of crosslink positions 
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y(s a ,i) = y a ,i (taken from Eq. ([2])), and gives 

W (&s a )= min W b [y a ]. (5) 

y(«a),y(»a,0=ya,< 

Technically, this is achieved by performing a cu- 
bic spline interpolation through the set of points 
{(sa,i,ya,i)}i=i,...,n a i- This can be shown to be equiv- 
alent to the minimization of the bending energy in 
Eq. (0J [3(J. As d 3 y ati /ds 3 is proportional to the trans- 
verse force in the fiber, the discontinuities of the cubic 
spline in its third derivative reflect the external trans- 
verse force that is needed to keep the fiber in its deformed 
shape. 

Now, we assume that neighbouring secondary fibers are 
free to react to the movement of the primary fiber. This 
allows the fiber deflection y a to be different from y«,i 
and may lead to a reduction of the bending energy on the 
primary fiber, however, at the cost of deforming the sur- 
rounding, i.e. by spreading the mode to the neighbouring 
filaments. The amplitudes of the so generated secondary 
floppy modes may be found from solving Eq. ([3]) for 6z ai . 
We thus find 



(6) 



which highlights the fact that a secondary mode of am- 
plitude Sz ai occurs when the actual transverse deflec- 
tion y a .i is different from the floppy-mode prescription 
y a .i- Furthermore, due to the scalar product with e a| ,-, 
the displacement y^ perpendicular to the plane defined 
by the two intersecting polymers does not contribute, 

Ya,i ■ g <M = 0- 

With Eq. ij^J) we find that Eq. ((5|) has to be modified 
by the bending energy contribution Wq from the neigh- 
bouring filaments giving 



Wi(& Q ) = min W b [y Q ] 
y(s=) \ 



E 

i=l 



W (Sz ai ) 



(7) 



Unlike in Eq. ([5|), where the crosslink variables on the 
primary fiber were constrained to be y(s a ,i) = y a ,ii here 
they remain unconstrained and move such that the to- 
tal energy, deriving from both primary and secondary 
fibers, is minimized. Note, however, that the deflec- 
tions on the secondary fibers are still constrained and 
given by Eq. ©, y{s auj ) = y a „j- This may be cor- 
rected for by taking into account further levels of fila- 
ments (tertiary, . . .), thus defining a sequence of energies 
(Wo, Wi,W2, ...Woo) the fixed point of which is found by 
substituting on both sides of Eq. ((7|) one and the same 
asymptotic function Woo- 

Since the resulting expression still depends on the 
quenched random network structure in a complicated 
way, we have recently proposed an effective medium ap- 
proximation that uses the averaged ([Woo]) = W in- 
stead [18j . For reasons that will immediately become 
apparent we have defined two averaging procedures. The 
angular brackets (.) denote averaging over the random 



variables on the primary fiber, the crosslink positions s a> i 
and angles 8 a ^. The probability distributions of these 
variables provide the most important characterization of 
the architecture of the network. The brackets [.] denote 
averaging with respect to the remaining randomness in 
the subsequent hierarchies of fibers. Mathematically, the 
effective medium approximation is implemented by inter- 
changing this latter average with the minimization oper- 
ation. Physically, this amounts to assuming that one and 
the same medium (W) is felt by all the crosslinks on the 
primary fiber. One thus arrives at the final equation 



W(Sz a ) = 



mm 

v y(*«) 



W b [y a ] 



E 



(8) 



where 6z ai is given by Eq. ©. In principle, Eq. |5| has to 
be solved self-consistently for the function W[x). Since 
we are concerned with small displacements only, the en- 
ergy may be expanded to harmonic order as W(x) — 
kx 2 /2, which gives 

W(<fe Qi ) = ifcsin 2 6> Qjl (y a>i - y a ^f , (9) 

where we defined y a .i = y a ,i ■ & a ,i and similar for y a ,i- 
With this parametrization Eq. ([5]) has to be solved for the 
single unknown parameter k. 

Eq. ([8]) can be interpreted as follows. The total energy 
stored in the network upon axially moving a single fiber 
the amount Sz a has two contributions. The first term, 
corresponding to the bending energy of the primary fiber, 
Wb, dominates if the crosslinks follow the local floppy- 
mode associated with Sz a , such that 5z ai ~ {y a ,% ~ 
Va,i)- On the other hand, the energy is mainly stored in 
the surrounding medium if the crosslinks deviate strongly 
from the floppy-mode, y a ^ « 0, in which case the bending 
energy vanishes, W b ~ 0. Since medium deformations 
can only occur in the form of floppy-modes, the stiffness 
k of the medium is the same as the stiffness of the fiber. 
This allows to solve the equation self-consistently, which 
can easily be done numerically as will be explained in 
the appendix. There, we will also solve Eq. ([5]) for some 
exemplary network structures. 

It is worth mentioning that Eqs. ([5]) and ^ may be 
interpreted as the zero temperature limit (or the saddle- 
point approximation) to a fluctuating stiff polymer in a 
random array of harmonic pinning sites with stiffnesses 
given by fcsin 2 # a)i (see Fig. Compared to the "bare 
floppion" defined by the Eqs. ([5]), the excitation given 
by Eq. ([5]) is "dressed" and incorporates the interactions 
with the medium on a Cayley-tree level. 



B. Elastic Modulus 

In principle, the clastic modulus can be found by mini- 
mizing the energy, consisting of contributions of the type 
of Eq. ([5]) from each of the Nr fibers, with respect to 
the variables 5z a [44j. Compared to the full problem of 
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having to minimize the energy with respect to all degrees 
of freedom, that is all N c \ ~ Nfn c \ crosslink coordinates, 
this is only a minimization with respect to Nt <C N c \ 
variables. Still, this poses a challenging quenched disor- 
der problem which can only be tackled numerically. 

Here we reduce the calculation to an effective single- 
fiber problem, by making a simplifying assumption about 
the magnitude of the individual 5z a . We assume, that 
the fiber centers-of-mass r" m = (X a ,Y a , Z a ) follow the 
macroscopic strain field in an affine way, just as the cen- 
ters of the hexagons did in the honeycomb lattice. This 
is equivalent to assuming that the displacement field is 
affine on the scale of the fiber length If. Note, how- 
ever, that this does by no means imply that the elastic 
elements themselves undergo affine deformations, as will 
become clear below [45j |. For a given macroscopic shear 
7 = we find <5r" m = jY a e x and thus 



8z a = jY a cos < 



(10) 



which is just the projection of the affine displacement on 
the fiber axis, oriented at an angle <fi a G [— 7r/2, 7t/2] with 
respect to the x-axis. Using Eqs. (fTO]) and ([2|) one can 
write Eq. (J3J) as 



y Q:i = -8z™ cot6> Qjl e Qii + y Q 



(11) 



where we have defined 



Ma =7 Pa cos <Pn ~ Y ai '- . (12) 

\ COS0 Q ,j J 

Upon comparison of Eq. (|11[) with Eq. @ one may inter- 
pret (fe™ 1 as specifying the movement of the primary fiber 
relative to its surrounding. Note, however, that this rel- 
ative displacement <fe^ cl depends on the orientations <fi a 
and 4> a + a .i of the primary and the secondary fibers, 
as well as on the arc-length along the primary fiber (via 
Y ai ). In contrast to Eq. (J2]), which follows from mov- 
ing the primary fiber in a fixed environment, Eqs. (jlip 
and (|12p are derived from a joint movement of all fibers. 
For the following we are only interested in the typical 
magnitude of <fe„ el , which may be obtained by averaging 
over the angles <f> a and estimating the typical distance 
between the center of masses of the intersecting fibers as 
Y a — Y ai ~ If. We thus find that &£ el cx jlf. 

By assuming affine displacements of the fiber centers, 
we have thus succeeded in reducing the many-body prob- 
lem of the movement 5z a of Nf interacting fibers to the 
case of a single fiber moving the amount 5z™ 1 ~ 7Z/ rela- 
tive to its surrounding. The modulus can thus be calcu- 
lated from the knowledge of the energy W(Sz Icl ) calcu- 
lated in the previous section. From the definition of the 
modulus we find that G~/ 2 /2 = NfW/V, where V is the 
volume of the system. 

It should be made clear that the assumption of affine 
displacements of the fiber centers is different from the 
usual approach of assigning affine deformations on the 
scale of the single polymer segment 0, [lH, [l3[ . The latter 
would lead to deformations <5 a g oc jl s , proportional to the 



1 



FIG. 6: The energy in Eq. ((8| is that of a stiff polymer at- 
tached to springs of variable stiffness. It may therefore be in- 
terpreted as a polymer in a random potential. The potential 
is attractive and localized at pinning sites given by Eq. ([2]). 



length I s of the segment. Instead, axial displacements of 
the fiber as a whole are, by construction of the floppy 
mode, directly translated into non-affine deformations 
<5na cx jlf, which do not depend on the length of the 
segment but rather on the scale of the fiber length I f . We 
would like to emphasize the subtle difference betweeen 
"affine displacements" of single points (the fiber centers- 
of-mass), and "affine deformations" of fiber segments of 
length l s . 



IV. RANDOM NETWORK IN 2D 

Having presented the general concepts of the floppy 
modes and their energy we now proceed to introduce a 
simple model system where the ideas may be tested. The 
random two-dimensional network, the "Mikado model" 
[lj|, has the advantage that it only needs one structural 
parameter, the density of fibers p. 

The network is defined by randomly placing N elas- 
tic fibers of length !j on a plane of area A = L 2 such 
that both position and orientation are uniformly dis- 
tributed [31], H2, HH • The fiber-fiber intersections are as- 
sumed to be perfectly rigid, but freely rotatable crosslinks 
that do not allow for relative sliding of the filaments. The 
randomness entails a distribution of angles 9 6 [0, tt] be- 
tween two intersecting filaments 



P(0) 



(13) 



while distances between neighbouring intersections, the 
segment lengths l s , follow an exponential distribution (34| 



P(h) = (Is) 



-i e -W<Q 



(14) 



The mean segment length (l s ) is inversely related to the 
line density p — Nlf/A as (l s ) — ir/2p. The segments 
are modeled as classical beams with cross-section ra- 
dius r and bending rigidity k. Loaded along their axis 
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( "stretching" ) such slender rods have a rather high stiff- 
ness (l s ) — Ak/I s t 2 , while they are much softer with re- 
spect to transverse deformations kj_(l s ) = 3k/1^ ("bend- 
ing"). 

Numerical simulations [H, E3 for the effective 
shear modulus G of this network have identified a cross- 
over scaling scenario characterized by a length scale 

Z = hWf)~" ( 15 ) 

with « 2.84 [l6| [46| that mediates the transition be- 
tween two drastically different elastic regimes. For a fiber 
radius r>^ the system is in an affine regime where the 
clastic response is mainly dominated by stretching defor- 
mations homogeneously distributed throughout the sam- 
ple. The modulus in this regime is simply proportional to 
the typical stretching stiffness, G a s oc k\\((l s }) and inde- 
pendent of the fiber length I f . This is in marked contrast 
to the second regime at r <C £■ There, only non-affinc 
bending deformations are excited and the modulus shows 
a strong dependence on fiber length, 

and thus on density, G na oc Sp^ where p, = 2v + 1 « 6.67. 

As this latter non-affine regime is characterized by a 
ratio ku({l s ))/k±({l s }) ~ ((l s )/r) 2 ^> 1, and therefore a 
bending mode that is soft as compared to the stretch- 
ing mode, we may apply the floppy-mode picture de- 
veloped in previous sections to calculate the exponent 
p. To this end, we numerically solve Eq. (JSJ) for vary- 
ing numbers n c \ ~ plf of crosslinks per fiber. The av- 
erage (.) is thereby defined in terms of the probability 
distributions of Eqs. (fT5|) and (fT4")) . As a result, we find 
the fiber energy to scale as W ~ n^n/lj and a; « 5.75 
(see Fig. [7]). The shear modulus is infered from IF as 
G = 2pW/lfj 2 ~ p 6 ' 75 , which reproduces the exponent 
/i as measured in the simulation to a remarkable accu- 
racy. 

Based on the formalism of the preceding sections we 
have also developed [l8| a scaling argument that al- 
lows approximate solution of Eq. (jSJ) in terms of a sin- 
gle length-scale l m i n , which on a microscopic scale gov- 
erns the coupling of the fiber to the matrix. Since the 
stiffness k± ~ n/l^ of the individual polymer segment 
is strongly increasing with decreasing its length l s , we 
assume that segments with Z s < / m ; n rather deform the 
surrounding medium than being deformed itself, while 
longer segments l s > Z m i n are not stiff enough to deform 
the medium. The scale l m [ n therefore plays the role of a 
minimal length below which segments are stiff enough to 
remain undeformed. 

In terms of the crosslink deflections j/i, this implies 
that long (and soft) segments have yi ~ yi, while short 
(and stiff) segments have crosslinks that are in their orig- 
inal position yj « 0. Since the energy of a segment of 
length l a can be written as w(l s ) ~ k±y 2 ~ nyf/1 3 ., we 
find that the elastic energy is reduced by the amount 




FIG. 7: Solution of Eqs. (|8]) and <(9j for various numbers 
n c i of crosslinks per filament. The randomness is defined by 
Eqs, (|13|l and (|14[) . The different symbols at given n c \ relate to 
ensembles of varying size N = 100 . . . 1000. The lines W ~ 
and W ~ n®i serve to illustrate the quality of the fit. 

w(l m in) ~ K Vi/lmin as compared to the situation where 
also the short segments are deformed. In turn, the en- 
ergy in the neighbouring fiber is increased, where a floppy 
mode of amplitude Sz ~ is excited. The length-scale 
^min can therefore be determined by equating the energy 
reduction in the small segments, w(l m i n ), with the energy 
increase due to the additional floppy mode in the neigh- 
bouring fiber. This latter contribution can be calculated 
as an average over all segments of length l s > ^ m j n thus 
giving 

dl s P(l s )w(l s ) = w(l min ) . (17) 

As a result, we find ^ m i n ~ 1/ p 2 lf and thus for the average 
fiber energy W ~ K(plf) e /lf. This corresponds to an 
exponent /i — 7, which confirms the previous analysis. 

From Eq. (|14|) one may also induce a probabilistic 
interpretation of the length-scale / m i n . Segments with 
lengths l s < l m in will occur on average only once along a 
given fiber. This may be seen from solving the equation 

/ dl s P(l B ) , (18) 

stating that small segments will occur once in every 
n c i ~ plf crosslinks. There will therefore be typically 
one segment per fiber in the undeformed configuration 
yi « 0, while all others follow the floppy mode. 

These scaling arguments also provide additional in- 
sights into a more microscopic understanding of the 
crossover from affine to non-affine elasticity in random 
rod networks [161 . Upon rewriting the crossover scale £ as 
£ ~ I f (lmin /I f) 3 ^ 2 the scaling variable x = r/£ of Eq. (2) 
in [lj| takes the alternative form x~ 2 ~ k\\(lf)/kj_(l m i n ) 
such that the crossover scaling law of the modulus reads 

G(r,p)=p» 9(^(1 /)/fex(Un)). (19) 
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The scaling argument now identifies a competition be- 
tween the relative stiffness of the stretching and bending 
modes on the scale of the whole polymer fiber as the driv- 
ing force of the afhne to non-affine crossover. For large 
scaling arguments, k»(lf) 3> k±(l m - ln ), bending is the 
weaker mode and yields a non-affine response in the form 
of floppy modes. Stretching deformations become dom- 
inant as soon as fey (7/) becomes smaller than k±(l m i n ); 
this happens if the rigidity scale ? m i n < (r 2 ^/) 1 / 3 . 

The only requirement for the presence of a bending 
dominated regime (beyond the scale separation ku/k± 3> 
1) is a low coordination number, which for the random 
fiber network can be calculated as z = 4(1 — (plf)^ 1 ). 
This places the network below the rigidity transition 
for any finite If, while increasing the filament length 
If — > oo the critical coordination of z c = 4 is asymp- 
totically reached. As an implication the bending mode 
must eventually be suppressed. 

The above analysis clearly shows that the proposed 
floppy mode concept can be utilized to understand the 
bending dominated elasticity in the random fibrous net- 
work. It allows to extract the length-scale Z m i n that is ul- 
timately responsible for the strong density-dependence of 
the elastic modulus as found in the simulations. Most im- 
portantly, the length-scale i m i n is a special feature of the 
random architecture studied here. Other network struc- 
tures will not necessarily feature the same length-scale 
even though the basic formalism of the floppy bending 
modes can still be applied. The exponents characterizing 
the elastic response will thus depend on network archi- 
tecture, a fact which is also exemplified in the Appendix. 

In Ref. [27| we have furthermore applied the theory to 
explain the mechanics of reconstituted actin networks, 
where filaments are crosslinked and bundled by fascin. 
By taking into account the fact that bundles have to 
be characterized by a length-dependent bending rigidity 
k(L) [H, [H, H3] it was possible to explain the observed 
dependence of the elastic modulus on actin and fascin 
concentration. 



A. Nonlinear elasticity arising from geometric 
effects 

Here, we report on additional simulations probing the 
nonlinear modulus of the structure. Note, that in these 
simulations the material properties of the fibers remain 
linear, such that the nonlinearities result from geomet- 
rical effects only. As one can infer from Fig. [8] the net- 
work is strongly stiffening already at very small values 
of strain. Similar results have recently been reported 
in [38(, where the stiffening behaviour was attributed to 
a crossover from bending to stretching dominated elastic- 
ity. The floppy mode picture allows to give this crossover 
a microscopic explanation. As argued in Section lllBl the 
floppy modes of the fibrous network are only adequate for 
infinitesimally small displacements Sz. The construction 
embodied in Eq. @ keeps segment lengths invariant to 
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FIG. 8: Nonlinear "modulus" cr/7 in the bending dominated 
regime (plf — 30) for various values of the aspect ratio a — 
r/lf. Inset: The stress increases linearly up to a strain of 
about 1%. Normal stresses quickly rise in magnitude and 
eventually are of the same order and proportional to the shear 
stress. 



first order in Sz only, such that any finite deformation 
will necessarily lead to stretching of the bonds. 

Note, that this stiffening mechanism is not mediated 
by non-linear material properties of the fibers but rather 
is of geometric origin and is due to the specific structural 
arrangement in the fibrous architecture. It is therefore of 
different nature than the stiffening mechanism inherent 
to single scmiflcxiblc polymers, where an applied tension 
can stretch the polymer only as far as there is stored 
length available [39(. 

In the nonlinear regime we have also measured the nor- 
mal stresses o~ xx and o~ vy that act perpendicular to the 
principal strain direction. We found (see inset to Fig. [8]) 
that these stresses can become of the order of the shear 
stresses a xy and have a negative sign indicating that the 
network "pulls in" during the course of the deformation. 
A similar effect has recently been observed in rheological 
measurements on F-actin networks (40j and rationalized 
in terms of the highly nonlinear entropic stretching re- 
sponse of single polymers. Note, that in our simulations 
the same effect occurs within a purely mechanical pic- 
ture, where no material non-linearities are present. It is 
explained with the fact, that the additional amount of 
contour length necessary to undergo a finite floppy mode 
can only come from pulling in the fiber ends. This is 
equivalent to a network contraction which leads to the 
observed large normal stresses. 



B. Nonstraight fibers 

In real networks fibers will never be perfectly straight. 
We have argued above that in this case the scale of the 
fiber-length If must be viewed as the length-scale over 
which the polymer remains straight. With this in mind 
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our theory also holds for networks where fibers are non- 
straight, as long as the undulation wavelength A ~ If is 
larger than the distance between crosslinks (l s ). 

In this section we investigate the effects of introduc- 
ing undulations with wavelengths comparable to the 
crosslink distance, A ~ (l s ). To this end we have manu- 
ally generated zig-zag fibers by randomly displacing the 
crosslinks by some maximal amount A • If . A similar anal- 
ysis has been performed in Ref. [38| , where a substantial 
decrease in the degree of non-affinity of the deformation 
field has been found. Similarly, we find that the system 
develops a new crossover to a regime of affine bending 
deformations (see Fig. [9]), where the modulus scales as 
G oc Sp 3 , a behaviour well known from bending domi- 
nated cellular foams [ill, H3, El • 




FIG. 9: Shear modulus G (inset) and scaling function g of 
Eq. (|20|l for various values of A. Collapse is achieved by plot- 
ting GA M//2 as a function x = SpA 1 ^ 2 and known exponent 
(i = 6.7 The asymptotic regimes show the scaling prop- 
erties of straight fibers, g(x) ~ x 6 ' 7 , and foams, g(x) ~ a; 3 , 
respectively. 

In this new regime the bending deformations come 
from pulling out the zig-zags similar to the pulling of 
thermally activated polymer undulations. We find that 
the curves may be scaled by using the same length-scale 
'miii ~ Sp that served as a lower cut-off in segment 
lengths. The modulus thus takes the following scaling 
form 

G(A>) = A^/ 2 5 (A/Z min ) , (20) 

where the scaling function has the limiting form g(x <C 
1) ~ x^/ 2 to eliminate the A-dependency. For large val- 
ues of the scaling variable x 3> 1 we have to recover 
the scaling properties characteristic for foams, giving 
g{x 1) ~ x 3 / 2 . This analysis highlights once more 
the fundamental role played by the length-scale l m - m in 
establishing the elastic response of the network. Here, it 
acts as a crossover scale, that mediates the transition to 
a foam- like bending regime at strong disorder A ^> i m i n . 

Note, that by introducing kinks in the contour of the 
fibers, the floppy modes start to spread beyond the single 



fiber to which they were confined originally. A kink is 
most conveniently characterized by the angle ?/> through 
which the direction of the fiber changes at the location 
of the kink. By displacing a crosslink by the amount 
A one thus finds for the angle shn/> = A/Z s , where l a 
is the length of the segment that ends at the crosslink. 
Exciting the fiber with a floppy mode of amplitude Sz, a 
finite kink-angle ip leads to the fraction Sz' ~ Sz sin ip ~ 
SzA/l s being coupled into the neighbouring fiber. At the 
crossover, defined by x = A/Z m ; n ~ 1, we therefore find 
that for a segment of length l s = l min , Sz'(l m - m ) « Sz. 
In this situation the floppy mode is transmitted to the 
neighbouring fiber without attenuation of its amplitude. 
Since segments of length /,„;„ statistically occur once per 
filament, the crossover point also marks the onset of a 
complete delocalization of the floppy modes. 



V. CONCLUSION 

We started our discussion with the assumption that 
the elasticity of stiff polymer networks is governed by the 
action of the bending mode. This assumption is based 
on the recognition that in systems where the persistence 
length is large, bending as compared to stretching is by 
far the softer mode. The respective spring constants are 
scale-separated and obey the relation fcii/fci ~ l p /l ^> 1. 

One immediate implication of this scenario is that 
polymer end-to-end distances have to stay constant, 
which necessitates deformations that are highly non- 
affine. We have characterized this non-affine deforma- 
tion field by constructing the floppy modes of the struc- 
ture [HI]. These are defined as the set of crosslink dis- 
placements that do not lead to any stretching of bonds. 
With this microscopic deformation field it is possible to 
calculate the macroscopic elastic moduli on the level of a 
self-consistent effective medium theory that incorporates 
fiber-medium interactions within a Cayley-tree approxi- 
mation. 

As a result the anomalous scaling properties of the 
linear shear modulus as determined by computer simula- 
tions of two-dimensional random networks [la . Il7j | are ex- 
plained. The exponents are found to be a consequence of 
the special architecture of the network that features two 
different length-scales. On the mesoscopic scale the fiber 
length If induces a non-affine deformation field, with seg- 
ment deformations (5 na following the macroscopic strain 7 
as <5 na ~ 7//, instead of as S a g ~ -fl s , which would result 
from an affine deformation field. Microscopically, a sec- 
ond length l m i n plays the role of a minimal length below 
which segments are stiff enough to remain undeformed. 

We would like to emphasize that the construction of 
the floppy modes only relies on the presence of the meso- 
scopic length If, which is applicable to a broad class of 
networks. In the particular case of a random rod network 
we have found that the anomalous scaling properties of 
the shear modulus, previously found in simulations, cru- 
cially depend on the presence of a second length-scale 
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^min, which is a special property of this random architec- 
ture. The exponents found for random rod networks are 
therefore not immediately applicable to other systems. 
Having established the general theoretical framework, it 
is nevertheless straightforward to calculate the exponents 
for other types of networks in two and three dimensions. 
Indeed, we have applied the theory to reconstituted actin 
networks crosslinked and bundled with fascin, and found 
that the calculated exponents are in good agreement with 
the experimental results [27| . 

Finally, we also conducted simulations probing the 
non-linear elasticity of the random fibrous network as 
well as modified the network structure by introducing 
kinks in the contour of the polymers. The results con- 
firm the governing role of the identified length-scales and 
firmly establish that the non-affine floppy mode picture 
captures the essential physics of stiff polymer networks 
similar in spirit to affine rubber elasticity for flexible 
polymer gels. In view of this conceptual analogy, the 
next step could be to assess the importance of crosslink 
fluctuations, which have been neglected here (as in clas- 
sical rubber elasticity). By greatly reducing the number 
of fluctuating degrees of freedom to one per fiber (namely 
5z), the theory developed here may very well provide a 
new starting point for the analysis of the statistical me- 
chanics of stiff polymer networks. 
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APPENDIX A: SOLUTION OF EQ. © FOR 
VARIOUS NETWORK STRUCTURES 

In this appendix we provide some technical details on 
how to solve Eq. (jHJ) for various network architectures in 
two spatial dimensions. We will measure lengths in units 
of the fiber length // and energies in units of n/lf, where 
k is the bending stiffness of the fiber. Assuming harmonic 
energies W(x) = kx 2 /2 we rewrite Eq. (8J symbolically 
as 

k=(f(k,n cl ;{z i ,6 i })), (Al) 
where the function / is defined by 

(A2) 

and we used Eq. @ to substitute ?/,; = — cot 9iSz. 

The network structure enters Eqs. (|Aip and (|A2|) via 
the variables {zi,9i}, which relate to the locations 



of the crosss-links on the backbone of the primary fiber 
as well as the angles 9i between primary and secondary 
fibers. The ensemble average (.) can then be defined by 
the probability distributions P({9i}) and P({h}), where 
segment lengths are given by U = z%+\ — Zi. 

To illustrate the importance of structural features on 
the elastic properties of the network we solve Eq. (|A2j) 
for two types of distributions, relating to random and 
regular structures, respectively. The random network 
is characterized by probability distributions as given in 
Eqs. (flU)) and (Ti"4"|) . The regular network has only one 
segment length Iq = lf/(n c \ — 1) and an angular distri- 
bution similar to Eq. (fl"3)) but restricted to the interval 

[^mini 7T $min] • 

For a given realization of the randomness the function 
/ is calculated by performing the minimization with re- 
spect to the contour y(z). This is achieved in two steps, 
where first the bending energy is minimized for a 

given set of values {y(zi)}. As explained in the main 
text, this is equivalent to a cubic spline interpolation. 
The second step consists of a minimization with respect 
to the remaining variables {y(zi)}. 

Finally solving Eq. ()A1[) the fiber stiffness k*(n c \) is 
determined as a function of the number n c \ of crosslinks 
per fiber. A graphical solution for n c i = 40 for various 
network structures is presented in Fig. 1101 The function 
(f(k)) is plotted as function of k. The sought after value 
k* is found at the point of intersection with the bisecting 
curve. 
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FIG. 10: Graphical solution of Eq. (|A2[I for a fiber with n c i = 
40 crosslinks imbedded in networks of varying architectures. 
{/(&)) i s plotted as a function of k for ensembles of varying 
sizes. The solution to Eq. (IA2[l is found by intersecting the 
curves with the bisecting line (dashed). 

The different curves for a given network structure cor- 
respond to ensembles of varying size. They seem to di- 
verge in the limit k — > oo. In fact, in this limit only the 
bending energy Wb contributes to Eq. (|A2[) and yi » f/j. 
This may make the averaging procedure ill defined, for 
example in the case of Eq. (fi"4"|) where the segment lengths 
l s can become arbitrarily small. The resulting segmental 
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bending energy Wf, ~ l~ 3 shows a divergence and does 
not have a well defined average value. 

As one can see from Fig. [10] the resulting fiber stiff- 
ness k* very sensitively depends on the randomness in 
the segment lengths, while crosslink angles only play a 
minor role. This is made particularly clear by comparing 
the random-segment and the regular-segment network in 



terms of the exponent x, which is defined by k* ~ nj. As 
stated in the main text the random network has x ~ 6, 
while for the regular network we find x w 4. We have 
shown above that the former result derives from the pres- 
ence of the length-scale Z m j n . In contrast, the latter is 
simply obtained by calculating the bending energy of n c i 
segments each of length l , Wb ~ Ucik/Iq ~ n^. 



[1] M. Rubinstein and R. H. Colby, Polymer Physics (Oxford 

University Press, 2003). 
[2] W. Kuhn, Kolloid Zeits. 76, 258 (1936), F. T. Wall, J. 

Chem. Phys. 10, 132 (1941), P. J. Flory and J. Rehner 

Jr., J. Chem. Phys 11, 512 (1943). 
[3] A. Bausch and K. Kroy, Nature Physics 2, 231 (2006). 
[4] B. Wagner, R. Tharmann, I. Haase, M. Fischer, and A. R. 

Bausch, PNAS 103, 13974 (2006). 
[5] R. Tharmann, M. M. A. E. Claessens, and A. R. Bausch, 

Phys. Rev. Lett. 98, 088103 (2007). 
[6] M. L. Gardel, J. H. Shin, F. C. MacKintosh, L. Mahade- 

van, P. A. Matsudaira, and D. A. Weitz, Science 304, 

1301 (2004). 

[7] Y. Tseng, B. W. Schafer, S. C. Almo, and D. Wirtz, J. 

Biol. Chem. 277, 25609 (2002). 
[8] K. Kroy and E. Frey, Phys. Rev. Lett. 77, 306 (1996). 
[9] F. C. MacKintosh, J. Kas, and P. A. Janmey, Phys. Rev. 

Lett. 75, 4425 (1995). 
[10] C. Heussinger and E. Frey, Phys. Rev. Lett. 96, 017802 

(2006) . 

[11] C. Heussinger and E. Frey, Phys. Rev. E 75, 011917 

(2007) . 

[12] J. H. Shin, M. L. Gardel, L. Mahadevan, P. Matsudaira, 

and D. A. Weitz, PNAS 101, 9636 (2004). 
[13] C. Storm, J. J. Pastore, F. C. MacKintosh, T. C. Luben- 

sky, and P. A. Janmey, Nature 435, 191 (2005). 
[14] L. J. Gibson and M. F. Ashby, Cellular Solids: Structure 

and Properties (Cambridge University Press, Cambridge, 

1999). 

[15] E. Frey, K. Kroy, J. Wilhelm, and E. Sackmann, in Dy- 
namical Networks in Physics and Biology, edited by G. 
Forgacs and D. Beysens (Springer, Berlin, 1998), Chap. 9. 

[16] J. Wilhelm and E. Frey, Phys. Rev. Lett. 91, 108103 
(2003). 

[17] D. A. Head, A. J. Levine, and F. C. MacKintosh, Phys. 
Rev. Lett. 91, 108102 (2003); Phys. Rev. E 68, 61907 
(2003). 

[18] C. Heussinger and E. Frey, Phys. Rev. Lett. 97, 105501 
(2006). 

[19] Rigidity Theory and Applications, edited by M. F. 

Thorpe and P. M. Duxbury (Kluwer Academic/Plenum 

Publishers, 1999). 
[20] J. C. Maxwell, Philos. Mag. 27, 27 (1864). 
[21] V. S. Deshpande, M. F. Ashby, and N. A. Fleck, Acta 

mater. 49, 1035 (2001). 
[22] C. R. Calladine, Int. J. Solids Struct. 14, 161 (1978). 
[23] S. Pellegrino and C. R. Calladine, Int. J. Solids Struct. 

22, 409 (1986). 
[24] S. Pellegrino, Int. J. Solids Struct. 30, 3025 (1993). 
[25] S. Alexander, Phys. Rep. 296, 65 (1998). 



[26] S. J. Winder and K. R. Ayscough, Journal of Cell Science 
118, 651 (2005). 

[27] O. Lieleg et. al., Phys. Rev. Lett. 99, 088102 (2007). 

[28] R. Cotter, K. Kroy, E. Frey, M. Barmann, and E. Sack- 
mann, Macromolecules 29, 30 (1996). 

[29] J.-P. Collet, H. Shuman, R. E. Ledger, S. Lee, and J. W. 
Weisel, PNAS 102, 9133 (2005). 

[30] J. Stoer and R. Bulirsch, Introduction to Numerical Anal- 
ysis (Springer, 2002). 

[31] J. A. Astrom, J. P. Makinen, M. J. Alava, and J. Timo- 
nen, Phys. Rev. E 61, 5550 (2000). 

[32] J. A. Astrom, S. Saarinen, K. Niskanen, and J. 
Kurkijarvi, J. Appl. Phys. 75, 2383 (1994). 

[33] M. Latva-Kokko and J. Timonen, Phys. Rev. E 64, 
066117 (2001). 

[34] O. Kallmes and H. Corte, Tappi 43, 737 (1960). 

[35] M. Bathe, C. Heussinger, M. Claessens, A. Bausch, and 
E. Frey, |q-bio.BM/06 07040, 2006. 

[36] C. Heussinger, M. Bathe, and E. Frey, Phys. Rev. Lett. 
99, 048101 (2007). 

[37] M. M. A. E. Claessens, M. Bathe, E. Frey and A. R. 
Bausch, Nat. Mater. 5, 748 (2006) 

[38] P. R. Onck, T. Koeman, T. van Dillen, and E. van der 
Giessen, Phys Rev Lett 95, 178102 (2005). 

[39] J. F. Marko and E. D. Siggia, Macromolecules 28, 8759 
(1995). 

[40] P. A. Janmey, M. E. McCormick, S. Rammensee, J. L. 
Leight, P. C. Georges, and F. C. MacKintosh, Nature 
Mat. 6, 48 (2007). 

[41] A. M. Kraynik and W. E. Warren, in Low density cel- 
lular plastics, edited by H. Hilyard and C. Cunning- 
ham (Kluwer Academic Publisher, Amsterdam, 1994), 
Chap. 7. 

[42] This exhausts all possible floppy modes in two spatial 
dimensions, while in three dimensions additional modes 
are found that relate to the deflection y x of individual 
cross- links. 

[43] The probability of intersection of two fibers with rel- 
ative angle 6 is proportional to the excluded volume 
A = sin 6^/2. 

[44] By minimizing the energy, instead of calculating the par- 
tition function associated with the variables Sz a , one ne- 
glects fluctuations of the crosslinks. 

[45] Similarly, the elements of the honeycomb unit cell do not 
undergo affine deformations [3], even though the unit 
cell as a whole does. 

[46] The density Sp is measured relative to the rigidity perco- 
lation transition occuring at p c lf = 6.7 



